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Abstract 

K^ , Wave packets in a system governed by a Hamiltonian with a generic nonlinear 

f^ . spectrum typically exhibit both full and fractional revivals. It is shown that the 

Q ' latter can be eliminated by inducing suitable geometric phases in the states, by 

^1/ ' varying the parameters in the Hamiltonian cyclically with a period T. Further, 

On \ with the introduction of this natural time step T, the occurrence of near revivals 






'i 



can be mapped onto that of Poincare recurrences in an irrational rotation map of 
the circle. The distinctive recurrence time statistics of the latter can thus serve as 



j^ ■ a clear signature of the dynamics of wave packet revivals. 
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1. INTRODUCTION 

There has been considerable interest in the dynamics of wave packets right from the 
early days of quantum mechanics. ^-' In recent years this subject has gained further im- 
petus from a variety of experiments in molecular systems, ^-* Rydberg wave packets, ^-' 
semiconductor quantum wells, *^-' etc., which require a detailed understanding of wave 
packet evolution - e.g., investigations in molecular physics (using femtosecond laser pulse 
techniques) based on the phenomena of full and fractional revivals ^-' of a vibrational wave 
packet. 

Consider the evolution of an initial state |'?/'(0)) of a system governed by a Hamiltonian 
H . In general, if |'?/'(0)) is not an eigenstate of H ^ the correlation or overlap function 

c{t) = imoMm' (1-1) 

will decrease from its initial value of unity as t increases. Under special circumstances, 
however, C{t) may return to its initial value at some particular instant of time, and a 
revival is said to occur. One might therefore expect that revivals would be facilitated 
if the system is prepared in an initial state \ip{0)) which is a superposition (a "wave 
packet") of stationary states of H, sharply peaked about some energy eigenvalue -E„„. 
(We consider, for simplicity, the one-dimensional case.) In the linear case of an equally- 
spaced spectrum, it is easy to show that revivals are a consequence of a simple periodicity. 
In general, however. En is a nonlinear function of n. On expanding En in a Taylor series 
about Ent), it turns out that the revival times depend on the coefficients of the linear 
and quadratic terms in this expansion. ^^ Cubic and higher order terms lead to (rare) 
"super-revivals" ^-' and can be neglected in most realistic situations, if the wave packet 
is peaked sufficiently sharply about En^- The quadratic term in the Taylor expansion 
leads to so-called fractional revivals that could occur between two full revivals: the initial 
wave packet evolves to a state that can be described as a collection of a small number 
of spatially distributed sub-packets, each of which closely reproduces the initial state ^-'. 
Experiments with Nal molecules display complex wave packet evolution that suggests 
fractional revivals of molecular wave packets. ^'^°-' 

Revivals and fractional revivals of a wave packet are, of course, a manifestation of 
the interference between the constituent basis states, each of which acquires a different 
phase during its temporal evolution. In general, the occurence of fractional revivals is 
contingent on delicate arithmetic properties of the numerical values of the parameters 
that appear in the Hamiltonian, e.g., the closeness of certain ratios of these parameters to 
rational numbers. The question that arises then is whether one can make the phenomenon 
of revivals more "robust" by suppressing the plethora of fractional or partial revivals in 



favour of clearly-signalled full (or nearly full) revivals, even though the spectrum is non- 
linear. We shall show that this is indeed feasible, by exploiting the possibility of inducing 
geometric (Berry) phases in the states over and above the dynamical ones acquired by 
unitary evolution. As is well known, such phases may occur if at least two parameters in 
the Hamiltonian are varied cyclically (with a period T, say) and adiabatically (the origi- 
nal, simplest setting for Berry phases ^^•'). This variation can be tailored so as to cancel 
out all fractional revivals in a generic nonlinear Hamiltonian. The cycling of parameters, 
however, means that it is only at instants of time separated by an interval T that the 
Hamiltonian returns to its original self. Therefore, with the introduction of this natural 
time step T into the problem, it is only meaningful to speak of possible revivals of a state 
at the instants T, 2T, ■ ■ -. It turns out that nearly full revivals ("near revivals") are still 
possible after this discretization of time. Moreover, the statistical distribution of these 
events is found to be precisely that of Poincare recurrences in the rotation map on a circle. 
The plan of this paper is as follows: In the next section, we give a brief review of wave 
packet revivals with particular emphasis on fractional revivals. In Section 3 we show that 
by inducing suitable geometric phases in the basis states, we can eliminate all fractional 
revivals. The formalism is also illustrated explicitly. Finally, in Section 4 we indicate 
briefly how the statistics of near revivals in the reduced problem can be mapped onto 
that of recurrences in an irrational rotation of the circle. 

2. FRACTIONAL REVIVALS: REVIEW 

Consider a system with a time-independent hermitian Hamiltonian iJ, with spectrum 
{En} and eigenstates {|0n)}- Let the system be prepared in an initial state |'?/'(0)) that 
is a superposition of the {|0n)}, sharply peaked about some riQ. As mentioned earlier, we 
expand En as 

En = En, + (n - no) < + (l/2)(^ - no)'< + ■ ■ ■ (2.1) 

As we wish to analyze only revivals and fractional revivals, we retain only terms up to 
the second order in Eq. (2.1) and shift n by hq for notational simplicity, to arrive at the 
quadratic form 

En = Co + Cin + C2n^ . (2.2) 

The coefficients Ci evidently depend on the parameters that occur in H. We shall assume 
that Ci, C2 > 0. (The modifications necessary in other cases are easily worked out.) In 

the |0„)-basis the time evolution operator U{t) = exp [—iHt/h] has the representation 

00 
f/(t) = Y. exp[-^ (Co + Cin + C2 n')t/h] |0„) (0„ | . (2.3) 



For a full revival (C(t) = 1) to occur at time t, U{t) must reduce to the unit operator 

(apart from a possible overall phase factor), i.e., {Cin + C2n^)t must be an integer multiple 

of 271% for every n in the summation. The following cases arise: 

(i) Ci 7^ 0, C2 = (equi-spaced or linear spectrum): Revivals of an initial state occur 

with a period T^ev = 2tt%/Ci. 

(ii) Ci = 0, C2 7^ 0: Revivals occur with a period Trev = 2'Kh/C2- 

(iii) Ci,C2 7^ 0, C1/C2 = a rational number r/s: Once again, full revivals occur with a 

fundamental revival time Tj-ev = 2Tihs/C2- A specific example is provided by the Hamil- 

tonian a) a^ = a'^a{a'^a — 1) that is relevant to wave packets propagating in a Kerr 

medium. ^^' It is evident that, in this case. En is proportional to n{n — 1) which is an even 

integer for every n. 

(iv) Ci,C2 7^ 0, C1/C2 irrational (the generic case): As the condition {C2n^ + Cin)t = 

2'n'hm {m = integer) cannot be satisfied for alln at any value oft, full revivals are no longer 

possible. However, at certain instants of time the quantity C{t) could come arbitrarily 

close to unity, producing a near revival. 

Between occurrences of full (or near) revivals, the wave packet breaks up into a finite 
sum of subsidiary packets at specific instants of time, provided the spectrum is nonlinear, 
i.e., C2 7^ 0. ^' These fractional revivals occur at times t given by t = 7rhr/C2 s where r 
and s are mutually prime integers. It can be shown that at these instants the evolution 
operator U can be expressed as a finite sum of operators Up, in each of which the phase 
factor multiplying the projection operator |0„) (0„| is linear in n: That is, 

Ui7vhr/C2s) = J2 4"'^ Up, (2.4) 



with 



and 



<'^^ 



p=0 



i-i 
'I /I) J2 exp [-(ircPr/s) + {2i7rkp/l)] , (2.5) 

k=0 



Up = J2 expi-tnOp) |0„) (0„| , Op = n [{Cir/C2s) + {2p/l)] . (2.6) 



n=0 



It is this decomposition of U which is responsible for fractional revivals, as can be seen by 
investigating the action of the operator Up on the initial state. For instance, if the initial 

state |'?/'(0)) is the "coherent state" 1^;) given by 



\z) = exp {-\z\y2) ^ |0„) , ^ G C , (2.7) 

n=0 V ^' 



Eqs. (2.4)-(2.6) yield 



i-i 
{nhr/C2s)) = exp (-ivrCor/Css) ^ 4"''M'2exp(-i^j,)). (2.8) 

p=0 



The state at time -v:h,r jC^s is therefore a weighted sum of wave packets, each of which has 
the same form as the initial state 1^;). 

It is clear from the foregoing that (i) fractional revivals arise from the quadratic 
dependence of En on n, and (ii) a whole host of such revivals of varying intensities can 
appear in a given case, depending on the numerics, i.e., the precise values of C^ and the 
integers r and s for which these revivals are detectable. In the next two sections we show 
that, by inducing suitable geometric phases in the basis states |0n), we can eliminate these 
partial revivals and restore near revivals at specific instants of time. 

3. SUPPRESSION OF FRACTIONAL REVIVALS 

We now consider the situation in which the Hamiltonian contains a set of "slow" parame- 
ters R that can be varied adiabatically and cyclically with period T, as in Berry's original 
setting for the geometric phase: i.e., Ry = Ro- 

In order to keep track of the parameter dependence in the basis states as well, we 
denote by |0„,R) the n eigenstate of the instantaneous Hamiltonian iJ(R). Then, at 
the end of each cycle of period T, |0„, R) picks up a geometric phase 7„ (over and above 
the i?„-dependent dynamical phase factor), given by the expression ^^^ 



7n 



^(0„,R|(Vr|0„,R)) ■ rfR, (3.1) 



where the integral runs over the corresponding closed loop in parameter space. Substi- 
tuting for |0„,R) at the end of a cycle, the time evolution operator at time T therefore 
becomes 



f/(T) = ^exp 



T'ln - (V^) / dtEn{Il. 

Jo 



Ro)(0„,Ro|. (3.2) 



To proceed, we need to know the n-dependence of 7„. For a general nonlinear spectrum 
En, one may expect a dependence of the form 

7„ = Oo + em + Gzn^ + ••• (3.3) 

This is generic, as it only requires that 7„ be a regular function of n (see Eqs. (3.13), 
(3.14) ff'. below). As in the expansion of En, only terms up to Oln^) in Eq. (3.3) are 
relevant for our present purposes. The coefficients ©« will clearly depend on the manner 
in which the parameters R are varied. We now substitute for 7„ and En in Eq. (3.2) 
from Eqs. (3.3) and (2.2), respectively, taking into account the fact that the coeffcients 
Ci are now time-dependent owing to the variation of parameters in H. Moreover, in each 



cycle of period T the state |0„, R) picks up the same additional geometric and dynamical 
phase. Hence the time evolution operator at time kT (where /c = 1, 2, ■ ■ ■) is given by 

U{kT) = ^exp iA;(z/o + z/in + z/2n^) |0„, Rq) (0„, Ro| (3.4) 

n 

where 

u, = Bi-il/h) r ait) dt,{t = 0,1,2). (3.5) 



/o 

From the discussion in Section 2, it is clear that all fractional revivals of a wave packet 
will be eliminated if the coefficient z/2 of n^ in the foregoing expression for U{kT) vanishes: 
this happens if we arrange the variation of the parameters such that 

62 = (l/h) f C2{t)dt. (3.6) 

Jo 

Once this is done, the exponent in U(kT) has only terms that are linear in n. Hence, 
at the relevant times kT {k = 1, 2, ■ ■ ■), the wave packet will no longer exhibit fractional 
revivals. 

To see in a little more detail how a geometric phase 7„ of the desired form may be 
generated, ^^^ consider first the classical one-freedom Hamiltonian 

H'{x',p') = Ap'^ + CV{x') (3.7) 

where A and C are positive constants, and V{x') is a potential that supports bounded, 
periodic motion. We assume that the quantum mechanical version, with a self-adjoint 
Hamiltonian H', has a non-degenerate, discrete spectrum {En} with normalized eigen- 
states (Ixn)}- The position-space wavefunction {x'\xn) = Xn{x') can be chosen to be a 
real function, from which it follows that the Berry phase 7„ that one may expect from 
the possible variation of A and C vanishes identically. Reverting to the classical case, 
consider now the canonical transformation {x',p') — > {x,p) where 

X = x' , p = p' - (B/A) f{x') (3.8) 

where i? is a constant and the function f{x) is to be specified. The transformed Hamil- 
tonian is 

H{x,p) = Ap" + B[pf{x) + f{x)p\ + CV{x) + {B''/A)f{x), (3.9) 

where we have written the cross terms in symmetric form in anticipation of quantization. 
Let f{x) be chosen such that the phase trajectories corresponding to H continue to 



represent bounded, periodic motion. Quantum mechanically, H is obtained by the action 
upon H' of the unitary operator 



W 



exp {-iBhF{x)/A) , F{x) = T f{x)dx. (3.10) 



Provided that this leads to a self-adjoint Hamiltonian H , it follows that H and H' are 
isospectral. The normalized eigenstates of H are given by |0„) = U\xn)- Under an 
adiabatic, cyclic variation of parameters R = (A, 5,(7), the Berry phase acquired by 

|0n) is 

7n = ^/(0n|(VR|0„)) ■ rfR. (3.11) 

As |0„) is normalized, this simplifies to 

In = I j{<pn\iyvJJ)\Xn) ■ dK. (3.12) 

Using the unitarity of [/, this yields 

In = {l/h)j {F)^Vn{B/A) -dR, (3.13) 

where 

(F)„ = / F{x)xl{x)dx. (3.14) 



Therefore the n-dependence of the Berry phase 7„ can be tailored by choosing the transfor- 
mation function f{x) appropriately, even though {-E„} remains the same for all acceptable 
choices of f\x). 

An explicit example is provided by the Poschl- Teller Hamiltonian ^^^ 

H' = Ap^ -Csech^x {A, C > 0) (3.15) 

which is unitarily transformed by W^ = exp[— {i B/ Ah) In cosh x] to H = W^ iJ' W^''" where 

H = Ap'^ + B[p{ta-nhx) + (tanha;)p] - \C+{B^/A)] sech^x + (B^A) . (3.16) 

Writing 77 = (1 + AC/Ah'^Y^'^, the spectrum of H is given by 

E„ = -Ah\r]-l~2n)y4: (3.17) 

where n = 0, 1, ■ ■ ■, [|(?7 — 1)] {[^] stands for the largest integer contained in ^). As En 
is quadratic in n, a wave packet whose evolution is governed by H will exhibit fractional 
revivals. Applying Eq. (2.4) to the case at hand, we can now identify the instants at 
which such revivals occur as nr/Ahs, where r and s are mutually prime integers. 



To induce a geometric phase 7„ in the eigenstate |0„) of if, it suffices in the present 
instance to vary just the two parameters A and B with a period T. It can be shown "'^^'' 
that the geometric phase in this case is precisely of the form 

7„ = 0o + em + 02^2. (3.18) 

Putting in the exact expression for ©2 in this instance, the condition (3.6) for the cancel- 
lation of fractional revivals then reads 

J[Vn{l/v^) X Vn{B/A)] ■ dS = h^ J^ A{t) dt , (3.19) 

where the integral on the left runs over the surface bounded by the closed loop in the 
{A, B) plane along which A and B are varied. 

Once fractional revivals are eliminated in the manner just described, we are left with 
linear n-dependences for both E^ and 7„. But, owing to the discretization of time in steps 
of T, this no longer implies trivial periodicity, i.e., that full revivals would automatically 
occur at integral multiples of T. However, as we shall show, near revivals will now occur 
generically without any further fine tuning of the parameters concerned. 

4. STATISTICS OF NEAR REVIVALS 

After the cancellation of the quadratic (n^) term in the phase factors in Eq. (3.4), the 
effective time development operator at time kT (/c = 1, 2, ■ ■ ■) is of the form 

U{kT) = exp(iA;z/o) ^exp(iA;nz/i) |0„,Ro) (0„,Ro| . (4.1) 

n 

Therefore an initial coherent state \ip{0)) = \z) given by Eq. (2.7) evolves to 

\^p{kT)) = \z exp(iA;z/i)) . (4.2) 

Correspondingly, the correlation function given by Eq. (1.1) becomes 

C{kT) = exp[2\z\'^{coskui-l)] . (4.3) 

Thus, if kui happens to be an integer multiple of 27r (recall that ui depends on T, cf. 
Eq. (3.5)), then full revivals occur with Trev = kT as the basic revival time. In general, 
however, ui is an irrational number (modulo 27r). Therefore, if 6q is the phase of the 
complex number z labelling the initial state |-?/'(0)), and 6k that of the corresponding 
number at time kT, Eq. (4.2) shows that the (discrete time) evolution of the state is 
entirely equivalent to the circle map 

Ok = 0k-i + 1^1 (mod 27r) (4.4) 



corresponding to a rigid "irrational" rotation. As is well known, ^^^ this map has no 
periodic orbits, is ergodic, and has a uniform invariant density. Given an e-neighborhood 
Je of the initial phase 6q, we have a Poincare recurrence to /^ at time kT if 6k E le- This 
corresponds to a near revival of the original wave packet, because it is easily shown that 
we then have 

C{kt) > l-\z\\^ . (4.5) 

The recurrence statistics for the rotation map is already known from certain "gap 
theorems" , ^^-^ and has found applications in the study of level spacings of oscillator 
systems ^^-' and recurrences in the coarse grained dynamics of quasiperiodic flow on a 
two-torus. ^^' In earlier work, '^^' we have discussed the application of these theorems to 
near revivals of coherent states of a deformed (displaced, squeezed) oscillator, and these 
results can now be taken over directly to the problem at hand. While the mean time 
between successive near revivals turns out to be simply 1-nT je in accordance with the 
ergodic theorem (the measure of /^ being tjl-n^ the probability distribution of this time 
interval is quite remarkable. It is typically concentrated at three specific values /ciT, k2T 
and (k\ + k'2)T ^ regardless of the values of v^ and e. The actual values of the integers k^ 
and A;2, and the relative frequencies of occurrence of the three different revival times, are 
of course dependent on v\ and e. ^''-' 

We have thus established an interesting link between revivals, anholonomies and recur- 
rences. These are manifestations, respectively, of quantum interference efffects, non-trivial 
topology in parameter space, and ergodic behavior in a discrete-time classical dynamical 
system. The distinctive nature of the generic distribution of near revival times that we 
have predicted should provide a clear signature, from the experimental point of view, for 
testing the validity of our formulation of certain aspects of wave packet evolution. 
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